Aug.27, 1998 



An Exact Universal Gravitational Lensing 

2 : Equation 

^— > . 
o . 

Simonetta Frittelli"' 6 and Ezra T. Newman 6 

^ '. a Dept. of Physics, 

i ■ Duquesne University, Pgh., PA 15282 

[>. ; & Dept. of Physics and Astronomy, 

3 ! University of Pittsburgh, Pgh., PA 15260 

O 

oo : 

Abstract 

o . 
cr 

^ ■ We first define what we mean by gravitational lensing equations in a 

bJO- general space-time. A set of exact relations are then derived that can be used 

>- ! as the gravitational lens equations in all physical situations. The caveat is 

that into these equations there must be inserted a function, a two-parameter 
family of solutions to the eikonal equation, not easily obtained, that codes 
all the relevant (conformal) space-time information for this lens equation 
construction. Knowledge of this two-parameter family of solutions replaces 
knowledge of the solutions to the null geodesic equations. 

The formalism is then applied to the Schwarzschild lensing problem. 
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I. Introduction 



The purpose of this note is to first define, on an arbitrary Lorentzian 
space-time, what we mean by exact gravitational lensing equations and then 
derive a version of these exact equations (which we believe should be of uni- 
versal applicability). Our definition reduces to the more usual lens equations^ 
when we consider approximations; namely small observation angle and per- 
turbation calculations off either flat or Robertson- Walker space-times. The 
basic idea will be to consider the one parameter family of past light-cones 
with apex on a time-like world-line that should be considered as the history 
of an observer; i.e., it defines what an observer can "see" and when the ob- 
server can see it. The gravitational lens equations are to be, loosely stated, 
the space-time positions of potential light sources, given in terms of the po- 
sition and time of the observation and the direction of observation (direction 
on the observers past tangent space light-cone). 

In Sec. II we will make precise what we mean by the "lens equation" 
and in Sec. Ill, using techniques obtained from V.I. Arnold'sBiJli theory of 
Lagrangian and Legendre maps, we derive a set of general relations, "lens 
equations", that with appropriate linearization reduce to the standard ap- 
proximate lens equation.El 

A caveat should be expressed here. First of all in any form of lensing 
equations one must be looking at the past null geodesies and one thus must 
have the metric tensor (needed only up to a conformal factor since only null 
geodesies are relevant) in order to calculate, from the geodesic equations, the 
null geodesies themselves. Though in our version, this information clearly 
must still be there, its specific appearance is circumvented by the assumed 
knowledge of a two-parameter family of solutions, u = Z(x a ,(X)i to the 
eikonal equation: 

g ab (x a )d a Zd b Z = 0. (1) 

The two parameters, ((, (), are the complex stereographic coordinates on the 
sphere. The (null) gradient of Z(x a ,(X) & t x a , namely p a = d a Z(x a ,(X)i 
sweeps out the sphere of null directions as (£, £) ranges over the sphere. 
The function Z(x a , (,(), which encodes all the conformal information of the 
space-time, is not easily obtained. It can be found exactly in spaces with 
sufficient symmetries, e.g., all conformally flat spaces, (Robertson- Walker), 
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Schwarzschild, Kerr, etc., but in most applications it will have to be found 
by perturbation methods. In any case, our lensing equations will be given 
explicitly and exactly in terms of Z(x a , (, () though when applied to physical 
situations, approximations will have to be made. 

Our starting point will be to assume that we have globally (perhaps given 
in different patches) the level surfaces of the function, u = Z(x a , (, (). 

Remark 1 For simplicity of discussion, we have assumed that the solutions 
to the eikonal equation will have been given by a function of the local co- 
ordinates, x a . Actually there might be singular regions (caustics) or self- 
intersections, where this is not possible; the surfaces could then be given in 
parametric form. We will nevertheless use the language of the simple func- 
tion, z(x a ,cX)- 

II. Defining the Lens Equation 

We first assume that our space-time allows the null geodesies to be in- 
finitely extendable into both the future and the past. In any region local 
coordinates will exist such that the world-line, £, of an observer can be given 
parametrically by 

x a = X a (r) (2) 
and the one parameter family of light-cones, £(Xq (r)), by 

x a = X a (XS(r), 77, fj, s) = X a (r, 77, 77, s) (3) 

where (77,77) label the sphere's worth of null directions at x a = Xq(t) and s 
is an affine parameter along each of the geodesies labeled by (r, 77,77). L a = 
d s X a 7^ is the (null) tangent vector to the geodesies. {Derivatives of X a 
with respect to the other parameters yield connecting vectors. They will 
not directly interest us here though they often are used in the study of the 
caustics of the light-cones, €(Xg(r)). We will approach the caustic issue 
from a slightly different point of view.} We further assume that the local 
coordinates x a are such that three of them, x l , are space-like and one, x°, is 
time-like. From this and from the non-vanishing null character of L a , using 
the implicit function theorem, one (at least one) of the space-like equations 
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(||) can be solved for s. Let us call that special coordinate x* and the other 
two space-like coordinates, x a . We then have that 

s = S(x*,T,i],fj) (4) 
and the remaining: three Eqs.© become, when s is replaced by S(x*, r, rj, rj) 

x° = X°(r,rj,fj, S(x*,T,rj,fj) = $°(x*,T,rj,fj), (5) 

x a = X*(T, V ,fj,S{x\T, V ,f}))=fr{x*,T,TI,f}). (6) 

If we think of a light source at the spatial position (x*, x a ) emitting light 
at time x°, then Eq.fl^) relates the time of emission x° at x*, with the time 
of arrival r and arrival direction, (77,77). We will refer to Eq.fl5|) as the time 
of arrival relation or function. 

The lens equation is defined by Eq.(|B|)/ it expresses two of the spatial 
coordinates, x a , in terms of the arrival time at the observer, r, the direction 
at the observer, that the observer sees the source and one of the source 
coordinates. Most often in applications the two coordinates x a are identified 
with some angular position coordinates (9, ip) of the source and the x* as 
some form of "radial" or "distance from observer coordinate", say D. In this 
case the lens equation reads 



9 = Q(,D,T,ri,f}), (7) 

Much of conventional lensing theory is based on finding approximate versions 
of Eq.(|7|). An important issue is the inversion of these equations, i.e., finding 
the observation angles in terms of the source position, namely when can one 
rewrite them as 

( V , fj) = (N(D, r, 9, if), N(D, r, 9, (p)) (8) 

and the multiplicity of these solutions. The condition for caustics on the past 
light-cone at D (the non-invertibility of (|7|)) is the vanishing of the Jacobian 
ofEq.(0); 
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J 



(9) 



We return to this issue in the next Section in a more general context. 

III. The Construction of the Past Light-Cone 

The construction of the light-cones, £(Xq (r)), Eq.(||) and hence, the time 
of arrival function and the lens equations, Eqs.@ and (|6]) was based on 
the integration of the null geodesic equations. In this section we will give 
an alternative construction of Eqs.@ and @ based on the solution to the 
eikonal equation. 

As discussed in the introduction, we assume that we have the two-parameter 
family of solutions, u = Z(x a ,(,Qi to the eikonal equation where for each 
value of (C ; C) the level surfaces of Z are null surfaces, i.e., d a Z is a null 
covector. Furthermore it is assumed that at each point x c , d a Z sweeps out 
the entire null cone at x a as (£, Q goes thru its range, S 2 . 

Remark 2 We repeat and emphasize that the level surfaces of the solu- 
tions to Eq.^) though referred to as "null or characteristic surfaces" are 
not strictly speaking surfaces; they can have self-intersections and in general 
are only piece-wise smooth. Though AmolS refers to them as "big -wave- 
fronts" we will continue to call them null surfaces. The intersection of a big 
wave front with a generic three surface yields a two-dimensional (small) wave 
front. 

The first thing that we want to show is that the light-cones, £(Xq(t)), 
with apex on an arbitrary curve, £, can be constructed from knowledge of 
the function Z. 

One sees immediately, from the eikonal equation, that the function 



for each fixed value of r,defines a two-parameter family of surfaces which all 
pass thru the point X£ (r) and which are all null surfaces. 



F(x a , r, c, C) = Z(x a , c, c) - Z(X%(t), c, C) = o, 



(10) 
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Remark 3 Note that d T F = -d T Z = -V a (r)d a Z ^ 0, with V a (r) = 
8 t Xq(t). The non-vanishing of d T F is because the scalar product of a non- 
vanishing null vector with a non-vanishing time-like vector is different from 
zero. It now follows from the implicit function theorem that Eq. (\ld{) is the 
implicit version of the function r = T(x a , (, (), i.e., we have 

z(x a , c, C) - TO, C, = o & r = T(x a , CO- (ii) 

The function T(x a ,(,C) ( or F) defines, what ArnolM% calls a generating 
family and is used to construct, via its envelope, the light-cone for each value 
ofr. 

Specifically the envelope of this family is constructed by demanding that! 

d ? F(a;V,C,C) = 0, (12) 
c\F(x a ,T,(,() = 

where (d^,d^ ) denote the derivatives with respect to the (CO- Assuming 
for the moment that fll2|) could be solved for the (C, C) = (r(^ a , t), T(x a , r)) 
then when they are substituted into ( |10| ) one obtains the function 

F(x a , t) = F(x a , r, T(x a , r),T(x a , r)) = Z(x a , T, T) - Z(X a (r), T, T) = 0. 

(13) 

Once again using the implicit function theorem with Eq. (|l2|) , Eq.(|l^ ), 
can be written as 

F(x a ,r) = Z(x a ,T,T) - Z(X a (r),r,f) = <£> r = f(x a ) (14) 

Using either F(x a , r) = or r = T(x a ) one can easily see that they define 
a one parameter, (r), family of null surfaces (and satisfy the eikonal equa- 
tion); this follows from either d a F(x a , r)| T = d a Z(x a , (, () or from d a T(x a ) oc 
d a Z(x a , C, C)- Furthermore, at x a = Xq(t), the gradients p a = d a Z(XQ (r), (, () 
sweep out all null directions. [At x a = Xq(t), Eq.(|12D is identically satis- 
fied and can not be solved for the (C, C) = (T(x a , r), T(x a , r)); all values of 
(Cj C) are allowed. See discussion below.] We thus see that F(x a ,r) = [or 
r = T(x a )} represents the family of null cones, €.(Xq(t)), with apex along £. 

The assumption that Eq. ([l2|) could be solved for (C, C) = (T(x a , r), T(x a , r)) 
depended on the non-vanishing of the determinant J of the matrix 
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Jij — 



d c d ( F c\d ( F 



c\c\F 



(15) 



J does vanish at the singularities of the "surface" F(x a ,r, 00, e -g-> & t the 



apex x a = Xq, among other regions. In fact the vanishing of J definesaa the 
caustics of the family of the cones, £(Xq (t)). 

In general, whether or not J = 0, Eqs.(|l^) and fllQ|) can be solved for 
other variables, namely some set of three (say x , which might be different 
in different regions) of the four x a , in terms of the fourth one (say x#), r and 
the (C,C), i-e-, 



x 



x I {x^ z 



(16) 



Note the important point that if the coordinates x a are such that three of 
them (say x J = {x a , x*}) are space-like and one of them is a time coordinate, 
a; , then Eq. (|T6|) has a stronger version, namely 



x° = ^(x*,r,C,0. (17) 
x a = ^(x*,r,CC) (18) 

where the two x J and the x* are the three space-like coordinates. That 
one can solve for the x° = X°(t, x*,(,() follows from the implicit function 
theorem and from the fact that T(x a ) satisfies the eikonal equation and hence 
df (x a )/dx° ^ 0. 

Eqs.(|I7D and QlJj) parametric representation of C(Xq(t)) via the 

null geodesies that rule it. For the different given values of the (CO, they 
are the null geodesies thru Xq. They are the explicit version of the implicit 
relations of Eqs.(|3) and (0). 

We thus see that Eqs.flTOD and ( p!2|) are equivalent to the lens and the 
time of arrival equations, © and (H) of the previous section. 

IV. Examples; Flat-Space and Schwarzschild Space- 
Time 

First, as a rather trivial example, using a family of Z(x a , (, () in flat space, 
we construct the family C(Xq (r)). A particularly useful two-parameter family 
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of solutions to the flat-space eikonal equation are the family of all plane-waves 
with can be represent eofi in the following fashion; 

Z(x a XX)=x a £ a (CX) (19) 

where d a Z = £ a {(,() represents the covariant version of the null vectors £ a , 
pointing in all possible directions, with Cartesian components given as 

f (C, C) = ^ ((i + CO, (C + OAl - 0, (CC - 1)) • (20) 

v2(l + CO 

Using the null basis set, (for each value of (C, ()), {£ a ((, C)> rn a ((, (), m a ((, (), 
n a ((, ()} where 



m a = (l + COV, m a = (l + C0<V a , (21) 



n 



a no, i / 1 i /•/•\2 



+ (i + C0 2 W a 



that have all vanishing scalar products except for £ a n a = —m a m a = 1. 
The lens and time of arrival equations, i.e., Eqs. (|T0|) and (pT2|), are then 



(x a -X a (r))£ a (C,0 = 0, (22) 
(x a -X a (r))m a = 0, (x a - X a (r))rn a = (23) 



where x a = X a (r), is an observers world-line, £. Eqs. (p2|) and (^), using 
(pT|), are easily solved for a parametric form of the light-cone as 



x a = X a (T) + s£ a ((,0, (24) 

with s a parameter along the null geodesic. If the observation is to be at an 
angle close to (say) the x* = x 3 — axis, then s = (x 3 — X 3 (t))/£ 3 ((, () and 
eliminating s from the remaining equations (|24|), we have for (|17|) and ([18]); 



x = X°(r) + (^ 3 -X 3 (r))^|l|, (25) 
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x 2 

the trivial time of arrival function and lens equation. A simple calculation 
of the Jacobian, 

J = d(x\x 2 )/d((,0 (27) 

verifies that it is different from zero and hence in the flat case there are no 
caustics. 

There are two ways to find the equations equivalent to Eqs. (p5|) and (f2~6|) 
in a Schwarzschild space-time; one could integrate the eikonal equation di- 
rectly, (this is easily doable by separation of variables!) or by integration 
of the null geodesic equations with appropriate initial conditions. (Actually 
this integration is most easily performed by the use of Hamilton- Jacobi tech- 
niques via the Eikonal equation.) In a recent work by Tom Kling, the exact 
Schwarzschild lensing equations were found.! We simply quote these results. 

Due to the spherical symmetry of the Schwarzschild lensing problem, with 
no loss of generality, one could place the observer at an arbitrary point on a 
fixed radial line - we chose the negative z axis - and place the sources in a 
special plane - the x — z plane. In polar coordinates that means the observer 
is at the angular position #0 = ir, O — 0. (For simplicity of discussion, 
since we working in a fixed plane, we let < 9 < 2ir, = 0, rather than 
< < 7r with = and ir.). The angular position of a source is then 
given by < 9 < 2tt and = 0. Rather then the usual Schwarzschild radial 
coordinate r, it is much more useful to use 

V2r 

and the retarded time coordinate, u = t — r — IMlogir — 2M) 

The time of arrival equation and the lens equation have the form 

u = Uo-U(l ,ip,l) (29) 
9 



X 1 (t) + (x 3 -X 3 {t)) 
A 2 (r) + (x 3 -A 3 (r)) 



i 1 (CO 
^ 3 (C,C) 
l 2 (C,C) 
^ 3 (C,C) 



(26) 



(28) 



e = e{ioM) (30) 

where (Z , Uq) are the position coordinates of the observer, who is assumed 
to be at the fixed position l on the negative z- axis. The coordinates of the 
source are (6, u, I) and ip is the observation angle, namely the angle between 
the incoming null geodesic projected into the rest-frame of the observer, 
and the inward directed radial vector at the observer, i.e., in the positive 
z direction. The roles of r and x*, in Eqs.(^) and ([Tj|), are played here 
by u and Z. The exact expressions for U and are given by the integral 
expressions, 



U = Uq - 



i di'd- 



sin 2 V> I' 2 f(l') ■ 



<o 2Pf(l')Jl 



sin 2 i{> I' 2 /(/') Jj 



W f(lo) 



i dl'(l- 



sin 2 ^ I' 2 f(l')~ 



2Pf(l') 



sin 2 V I' 2 /(/') 

P 

'(i 



(31) 



Wf(h) 



, lS sin-0 f l dl' smib f l dl' 



k Jio ^(i - sin y;;; o f ,} o *o >/7(w ^ */< 

where 



, _ sin 2 V> l' 2 f(l') * 

i? /Co) - 

(32) 



/m = i _ ^ = (i _ 2V2MZ) (33) 

r 

and Z is the positive root of the cubic equation 

(sin^) 2 Z' 2 (l - 2V2MI) -l 2 (l- 2V2Mlo) = 0. (34) 

If the source is situated on the positive z- axis, i.e., at 9 = then 
Eq.(|32f) becomes 



sin-0 dl' sinib f l dl 

= 7T rZ= / = ^=4 ' 



V v '0 /Go) 7 V V 'o /('o) 7 

(35) 
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which can be considered as an implicit function, F(ip, I, Iq)) = 0, that defines 
the observation angle if> when the source and observer are co-linear with the 
origin; if) is then the observation angle of the Einstein ring. 

Note that in this case we have cheated slightly in that we have used 
null coordinates for the time of arrival and lens equations rather than the 
spatial/time coordinates of Sees. II and III. 

Simple approximations! to Eqs. (|3l"l) , and (|35|) yield the standard 
linearized Schwarzschild time of arrival/lens equation and Einstein angle. 

V. Discussion 

In Sec. II, we gave a general definition of an exact time of arrival/lens 
equation that, in principle, is applicable to all physical situations. Its ap- 
plication depended on knowing how to construct the past light-cones from 
an observers world-line or on the construction of a two-parameter family of 
solutions to the eikonal equation. 

We are studying the possibility of applying these techniques to statistical 
perturbations of the homogeneous, isotropic cosmological models for which 
the eikonal solutions are known. 

From the point of view of contact with observation, it appears extremely 
unlikely that the exact equations discussed here, and their associated correc- 
tions from the linearized approach to lensing, implicit in the exact equations, 
will have observational consequences in the near future. Nevertheless, the ex- 
traordinarily rapid advances in observational techniques must make one du- 
bious of absolute statements that something is unobservable. To emphasize 
this point we end with a quote from Einstein's 1936 papei0 on lensing; 

• "....there is no great chance of observing this phenomenon." 
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